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Ssrstesmtlc  methods  of  obtaining  approximate  eolations  to  com¬ 
plicated  son- Umar  problems  art  the  subject  of  perturbation  theory. 
This  theory  probably  reaches  its  highest  development  in  celestial 
mechanics  vh ere  detailed  oowqy laxities  east  he  tikeg  into  account  in 
order  to  achieve  the  desired  accuracies.  Within  reoeat  years  there 
has  been  a  considerable  revival  of  interest  in  the  application  6f 
classical  Mechanics  theory  to  practical  problem  involving  satellites 
end  space  flight,  using  the  lasailton- Jacobi  theory,  ceaenloal  vari* 
ables,  and  von  Zelpel's  perturbation  Method  to  provide  seme  very  el* 
agent  solutions  to  artificial  satellite  problem. 

However,  may  engineers  without  considerable  training  or  ex* 
parlance  in  the  subject  here  difficulty  in  following  this  work.  Prob- 
•  ably  eost  of  the  difficulty  arises  because  of  their  unfsailiarity  with 
both  the  mthasMitlcal  and  astronomical  Jargon  involved  (e.g.,  "deter- 
nlnlng  functions, N  "secular,"  "periodic,”  "canonic,"  etc.)  and  be* 
cause  of  the  complexity  introduced  when  several  degrees  of  free  don 
are  Involved. 

She:  wen  Zslpel  perturbation  technique  at  the  tine  of  this  writing 
is  a  specialised  tool  available  principally  to  astronomers  and  those 
trained  in  celestial  aechaalca.  It  is  a  powerful  method  which  Should 

have  vide  application  to  problems  involved  In  other  fields-.  Accord*- * , 

_ _ — *' 

lngly,  this  Memorandum  presents  the  subject  from  an  engineering  rather 
than  mathaeatlcel  point  of  view,  using  relatively  simple  examples  to 
illustrate  the  mthod. 
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LIST  OF  SYMBOLS 


<0  or  k 

V 

E 

L 
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m 

PJ 
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angular  frequency  of  a  simple  harmonic  ays  tea 

generalized  coordinate  J  ■  1,  2,  3,  •  •  • 

the  Haailtoaian— -e^ual  to  the  am  of  the  kinetic  and 
potential  energies  of  the  system  (acmstimea  designated 
F  by  astronomers) 

canonical  variable  representing  momentum  conjugate  to 
the  angle  variable  t 

canonical  variable  representing  a  coordinate  angle  con¬ 
jugate  to  the  momentum  variable  L 

mass  of  a  particle 

generalized  momentum  J  -  1,  2,  3,  •  •  • 

a  email  positive  number  (6  <  1) 


Dote  above  symbols  represent  derivatives  with  respect  to  time 
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i.  ptnoDuenoM 

Itor  a»n  than  a  century  estrennm ere  and  nth— aticlaas  ban  u—d 
the  Emil  ton-  J aaobl  theory  and  related  perturbation  technique*  in  ex¬ 
plaining  and  predicting  the  notion  of  celestial  todies.  Many  of  the 
powerful  auxiliary  perturbation  methods  and  techniques  developed  by 
n—  such  as  Belaunay,  Kill,  Brown,  Shook,  and  iw  Zsipel  were  evelved 
in  solving  highly  specialised  problems  In  oelestlal  mechanics  related 
to  lunar  theory,  planetary  theory,  — d  eelar  system  stability.  These 
problems  vers  very  conplicated,  involving  three  (sad  sonatinas  four) 
degrees  of  freedom  and  entailing  an  — emeus  —oust  of  detailed  der¬ 
ivation  and  calculation.  Accordingly,  it  is  not  surprising  that  a*— e 
of  these  mathematical  tools  have  be—  familiar  over  the  yean  ealy  to 

t-  — 

persons  astronomically  oriented  and  indoctrinated  In  the  theory  of 
classical  mechanics.  Bsoantly,  wider  Interest  has  developed  as  en¬ 
gineers  and  applied  mathematicians  have  encountered  practical  problems 
in  space  flight  guidance. 

In  the  past,  it  was  the  astronomer  vho  was  Interested  In  pre¬ 
dicting  the  motion  of  some  celestial  body  for  an  extensive  period  of 
time,  say  50  or  100  years.  Today,  space  flight  or  satellite  systems 
are  being  planned  to  operate  tor  periods  ef  several  years.  The  pro¬ 
cedure  of  using  a  crude  mathematical  model  and  many  iterations 

on  an  electronic  computer  usually  breaks  down  in  problems  of  this  typo 
because  of  truncation  and  round-off  errors.  There  has,  therefore,  been 
a  need  fbr  the  vise  of  more  sophisticated  mathematics.  Within  the  past 
four  years  notable  contributions  along  these  lines  have  be—  made  by 
Brouwer,^  Garflnkel,^  Hori,^  Xozai,^)  and  many  others  using  the 
Hamilton- Jacobi  theory  and  a  perturbation  technique  usually  ascribed  to 
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von  Zeipel,^  although  others^' ^  contributed  methods  along  very 
similar  Hams.  Brouwer,  Hori,  and.  Oarfinkel  developed  s  complete  first- 
order*  solution  for  the  perturbed  motion  of  sn  artificial  satellite  with¬ 
out  air  drag.  Chaplets  second-order  eolations  of  this  problem  have  hoam 
developed  by  others.^  The*#  publications  appeared  in  l$9Pf  ****•  than 
a  of  papers  have  appeared  describing  the  adaptation  ef  ths  tadoipn 

to  various  perturbation  problems  in  celestial  mechanics  (see  references), 
ttafsrtunately,  many  of  thoao  psoblems  are  ecaplicated  in  dot  ail,  and  by 
the  very  nature  of  their  solution  rupiira  a  tedious  mount  sf  analytic 
derivation  and  development,  whereas  the  method  itself  la  siaplo  end 
frtsralgxtlbrvard.  A  thoxeutfi/  taBMledps  of  classical  mechanics  and  the 
calculus  of  variations  is  necessary  to  understand  the  underlying  theory 
involved,  but  only  an  elementary  acquaintance  la  necessary  in  order  to 
apply  the  technique  successfully  to  :wmgr  practical  engineering  problems. 

me  intent  hex*  Is  to  entlins  the  method  fw*  se  angtl— rfng^  vst*> 
or  than  a  mathematical,  point  of  view,  using  relatively  siaplo  examples 
involving  only  one  degree  of  fresdem  (l.o.,  the  harmonic  oscillator  and 
the  Van  der  Pol  equations) .  The  mathematical  theory  and  derivation  may 
be  found  in  texts  on  classical  or  celestial  mechanics. 

Because  of  the  similarity  of  the  von  Zeipel  perturbation  method  to 
that  of  Kryloff  and  Bcgpliuboff  method^12*1^)  vhich  is  well  known  in  the 
fields  of  physics  and  engineering,  oam  aectios  is  devoted  to  a  oonpariaoa 
of  the  two. 


•Referring  to  the  first  poser;,  of  the  coefficient  of  the  eeoond 
harmonic  of  the  earth** ^|«tsnUH~fbnetlea,  variously  designated  J2> 
k2,  B2,  etc. 
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II.  CABOHICAL  VAKEABLBB  AMD  TRAKSTOPMATIOKS 

<Sh«  von  Zeipel  method  consists  of  making  successive  mathematical 
trams  formations  of  variables  of  a  canonical  system  of  differential 
equations.  Ths  transformations  axe  performed  in  a  methodical  way 
'according  to  established  general  rules ,  so  that  the  final  solution 
is  obtained  in  a  certain  desired  fern.  For  example,  in  celestial 
mechanics  it  is  expedient  to  describe  the  variation  of  an  orbital 
parameter,  say  ,  In  the  form: 

-  qQ  +  ^(t-t^  ♦  long-period  trigonometric  terns 
+  ahort-perlod  trigpncsaetric  terms 

vhere  q^,  q^,  and  tQ  are  constants  of  the  motion.  B»e  linear  time 
varying  term  is  referred  to  as  being  "secular. "  B*e  long-period 
terms  may  involve  periodicities  very  great  in  magnitude  relative  to 
the  unperturbed  two- body  period,  whereas  the  short-ported  terms  are 
usually  com&ensurable  with  this  period.  To  an  astronomer,  it  is 
important  to  be  able  to  separate  mathematically  the  long-period 
from  the  secular  terms.  Iter  if  the  "long  period"  «•»  _««a>  •*] 

1000  years,  the  periodic  effect  would  be  difficult  to  separate  from 
the  secular  effect  by  observational  means. 

A  system  of  differential  equations  describing  a  mechanical 
system  Is  said  to  be  canonical  if: 


« 


Sh 


-Sh 


J  —  !• • >n 


where  H  is  the  Hamiltonian  equal  to  the  sum  of  the  kinetic  and 
potential  energies  of  the  system;  Pj  is  the  generalized  momentum 


(1) 
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and  q^  the  generalized  coordinate*  The  system  has  n  degrees  of  freedom. 
We  will  as  suae  that  B  may  be  written  so  that  it  does  not  contain  the 
time  explicitly.  That  is, 

H  ■  H(Pj,q^)  J  *  l...n 

Equations  (1)  are  a  set  of  2a  first-order  differential  equations.  Their 
solution  gives  and  q^  as  functions  of  the  tine  and  2n  arbitrary  con¬ 
stants.  As  an  example,  consider  the  simple  linear  spiring  where  a  par¬ 
ticle  of  mass  m  is  constrained  to  move  along  a  path,  its  displacmaent 

2 

frem  a  fixed  point  being  q.  The  spring  constant  is  k  . 

1  .2 

The  kinetic  energy  is:  ^ 

12  2 

The  potential  energy  is:  -g  k 


Designating  the  momentum  as  ■  mq^,  the  Hamiltonian  is: 

2 

Ho<px'«*>  "  s  +  i  4  (2) 

Hote  that  px  and  qx  are  both  functions  of  time  but  t  is  not  explicitly 
contained  in  Hq.  The  2n  canonical  equations  are 


dH 

•  o 

^  ’  3T 

X 


m 


(3) 


These  two  first-order  equations  may  be  readily  solved  for  q^t)  and  px(t) 
in  terns  of  2n  -  2  arbitrary  constants.  The  first  appears  to  be  rather 
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trivial,  merely  restating  the  definition  of  oonentun;  the  second  e (junc¬ 
tion  can  he  combined  vith  the  first  to  form  one  single  second-order 
equation  which  is  that  of  the  spring  or  simple  harmonic  oscillator. 


aq^  +  k2  ^  -  0  (h) 

Us  shall  he  concerned  hers  with  obtaining  solutions  to  the  first  order 
Xqs.  (3)  rather  than  Bq.  (>).  Ibr  the  perturbed  ease  (i.e.,  nose- 
linear  spring)  to  he  considered  in  the  next  section,  the  Haailtonian 
of  Bq.  (3)  is  not  in  a  form  amenable  for  perturbation  techniques.  A 
fundamental  theorem  in  the  transformation  theory  of  classical  mechanics 
is  that  a  transformation  of  variables  may  be  made  such  that  the  func¬ 
tional  form  of  the  Haailtonian  is  changed,  although  its  value  is  not, 
providing  t  is  not  explicitly  present  in  H.  If  an  appropriate  trans¬ 
formation  method  is  used,  the  differential  equations  describing  the 
system  will  roaain  canonical  although  of  a  different  functional  form. 
Bor  example,  let  us  assume  that  the  Hamiltonian  HQ  is  to  be  trans¬ 
formed  to  one  having  a  new  functional  form  with  new  variables  p^ 
and  q^.  Then 

Ho<VSc>  "  H>X>  (5) 


Sc 


(6) 


In  stationary  state  Hamilton- Jacobi  (B-J)  perturbation  theory,  the 
first  step  is  to  transform  the  Hamiltonian  of  the  unperturbed  prob¬ 
lem  (i.e.,  H  )  so  that  it  contains  the  "action"  variables  only  (i.e., 
o 

momenta,  1a  our  case).  In  practice,  the  matheaatics  of  this  step 
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can  usually  be  circumvented  by  Intuitive  reasoning  or,  better  yet,  by 
starting  with  ecmeone  else's  results  that  will  guarantee  the  desired 
functional  form  of  H  as  veil  as  the  B-J  equations  remaining  canonical. 
This  procedure  vlll  be  illustrated  later  in  Section  Y,  dealing  with 
applications  to  celestial  mechanics. 

Returning  to  the  exmpla  of  the  linear  spring,  the  Hamiltonian 
may  be  transformed  so  as  to  contain  only  a  momexrtw  variable  by  intu¬ 
itive  reasoning.  That  is 

HC(V  qx)_  H'(p^,  -  )  where  Hq  -  Hc'  (T) 

First  consider  the  solution  of  Eq.  (3)  with  the  conditions  m  •  1, 

%u>  "  0 

^j(t)  -  Qjo  cos  k(t-tD) 

(8) 

Px(t)  -  qLJC<t)  -  -  kq^  sin  k(t-tQ) 


At  time  tQ,  the  kinetic  energy  is  zero  and  the  total  energy  of  the 
system  is  equal  to  the  initial  potential  energy  or 


Ho- 


«*D- 


Ho 


•  total  system  energy 


(9) 


For  convenience,  ve  shall  designate  p^  end  its  conjugate  as  L  and 
l  respectively.  The  canonical  conditions  of  Eq.  (6)  and  the  desired 
transformation  of  Bq.  (7)  may  be  determined  intuitively  by  letting: 

k  v.2 

L  m 


2 


(ix>) 
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so  that 


l  m  k(t-tQ) 


H^P^)  -  H^(L,  -  )  -  k  L 


Then,  the  canonical  aquations  are 


or  L  -  const. 


K 


(11) 

(12) 


(13) 

(1*) 


The  solution  to  the  linear,  i.e.,  unperturbed,  problem  may  be  con¬ 
sidered.  as  auxiliary  equations  connecting  the  transformation.  That 
is 


cos  t 


(15) 


q^  »  -  /£  k"l  sin  l 


In  the  unperturbed  case,  L  and  the  angular  variable  l  are  "elements," 
or  "parameters,"  specifying  and  defining  the  simple  harmonic  motion 
of  Sq.  (15)  •  It  is  logical  to  assume  that  if  this  slsple  motion  is 
suddenly  perturbed  by  a  small  disturbing  force,  *qs.  (13)  could  still 
be  used  to  describe  or  approximate  the  resulting  motion]  however,  tbs 
elements  would  require  a  slight  adjustment  to  represent  the  motion 
accurately  for  a  given  short  Interval  of  time.  This  subject  is  dis¬ 
cussed  further  in  the  next  section. 
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III.  THE  VARIATIQH  OT  PARAMETERS  TECHNIQUE 

In  order  to  clarify  the  Material  vhleh  la  to  be  presented  in  the 
next  section,  it  may  be  worthwhile  to  review  briefly  the  variation  of 
parameters  technique.  In  the  preceding  section  the  solutions  of  Eqs. 
(3)  were  of  the  fora 

q^t)  ■  A  cos  (kt  +  B) 

<*> 

Px(t)  *  ^(t)  ■  -  A  k  sin  (kt  +  B) 

where  A  and  B  are  constants  of  integration  and  k  is  the  angular  fre¬ 
quency  or  in  terns  of  other  variables  which  have  been  used 

A  * 

B  «  -  kt 

o 

so  that  Sq.  (16)  corresponds  to  Eqs.  (8)  or  (15).  It  was  suggested 
that,  if  the  simple  harmonic  notion  described  by  these  equations  wee 
perturbed  by  a  small  disturbing  force,  these  solutions  could  still  be 
used  to  describe  or  approximate  the  resulting  notion.  However,  the 
elements  or  constants  A  and  B  would  require  slight  adjustments  to 
represent  the  motion  accurately  for  a  given  short  interval  of  time. 
Moreover,  as  time  progresses,  one  would,  expect  that  adjustments  to 
the  elements  would  be  continuously  or  periodically  necessary  to  achieve 
this  purpose.  Therefore,  these  "constants"  of  the  above  equations 
may  be  thought  of  in  the  perturbed  case  as  slowly  tine-varying.  In 
perturbation  theory  this  concept  is  usually  referred  to  as  the 


"variation  of  parameters"  method.  Generally ,  one  may  expect  that  the 
adjustments  to  the  elements  will  consist  of  additive  terms:  a  linear 
tine- varying  texm  plus  periodic  tens  of  Beveral  different  frequencies 
as  in  a  Fourier  series.  To  illustrate  this  technique,  consider  an 
equation  of  the  font 

‘4  +  V  t>  “  0  (17) 

where  the  solution  is  nearly  sinusoidal  or  quasi- harmonic  so  that 

Vx>  t)  -  k2  q^  ♦  l  X(qx,  px,  t)  (IB) 

In  the  shove,  t  is  a  parameter  characterizing  the  smallness  of 

2 

the  deviation  of  F  from  k  q^.  Bert,  if  ve  examine  the  cue  vhare  t 
is  not  explicitly  present  in  X,  Xq.  (17)  becomes 

4  +  *2  Sc +  ?  x<  v  px)  ■ 0  (19) 

where,  for  &  ■  0,  Eqs.  (l 6)  are  the  solution.  For  C  small,  end  positive, 
A  and  B  may  be  considered  nev  unknown  functions  of  time  and  determined 
in  such  a  manner  that  Eqs.  (16)  are  forced  to  be  solutions  of  Eq.  (19). 
To  do  this  ve  assume 

-  ^(t.  A,  B) 

where 

A  -  A(t)  and  B  -  B(t) 

so  that  differentiating  the  first  of  Eqs.  (16)  vlth  respect  to  time 


low  egcwinlng  the  above  equation  and  considering  the  requirement  for 
Bx(t)  ©f  Sqe.  (16),  it  is  necessary  that 

da 

»x(t)  -  -  -  Ak  sin  (kt  +  B) 

and 


»•  (kt  >B)  -  A§  sin  (kt  ♦  B)  -  0  (20) 

Similarly,  differentiating  px  of  Bqs.  (1 6)  with  respect  to  time  re¬ 
sults  in 

Px  *  \  -  -  k  sin  (kt  ♦  B)  -  Ak2  cos  (kt  ♦  B) 

-  Ak  ^  oos  (kt  +  B) 

ffieaminlng  the  above  equation  and  considering  the  requirements  for 
of  Zq.  (19)  and  for  (see  Sqs.  16),  it  is  necessary  that 

k  ^  sin  (kt  ♦  B)  +  Ak  f§  eos  (kt  +  B)  -  ftx(qjc,  px)  (21) 

Solving  Bqs.  (20)  and  (21)  'for  and  || 


^  sin  (kt  +  B) 
H  «  H  cos  (kt  ♦  B) 


(22) 


where 


x  »  x(  v  V 


11 


Bote  that  and  ^  being  proportional  to  t  vlll  be  slowly  time- 
varying  during  the  period  2*/k. 

Surnoarizing,  if  Eqs.  (22)  can  be  Integrated  so  that  A  and  B  are 
known  or  approximated  as  functions  of  time,  then  Eqs.  (l6)  would  rep* 
resent  a  solution  to  Eq.  (19)  •  In  other  words,  Eqs.  (l6)  can  serve 
to  represent  the  perturbed  t  t  o  solution  if  the  parameters  A  and  B 
are  considered  to  be  time- varying  according  to  Eqs.  (22).  Bow  re¬ 
turning  to  the  canonical  variables  L  and  1,  It  Is  clear  that  since 


then 


and  since 


then 


dl 

IF 


dA 

if 


B  ■  -  kt 

o 

i  -  k(t-tQ)  -  kt  ♦  B 


d  l 
3F 


•  k  + 


dB 

dF 


Bote  that  Eqs.  (22)  could  be  substituted  In  the  above  expressions  for 

Sr  jO  a 

and  ^  .  Therefore  to  obtain  a  solution  to  the  canonic  equations 

dL  5h 

as "  •  si 

dA  dB 

dF  “  3E 


where  E  is  the  Hssn.il tonian  of  the  non-linear  or  perturbed  problem,  is 
equivalent  to  solving  Eqs.  (22).  In  the  next  section  we  shall  con- 


12 

aider  von  Zelpel's  method  for  approximating  the  solution  of  the  aborre 
canonic  equations.  The  results,  of  course,  with  Eq.  (16)  (equivalent 
to  Eqs.  (8)  or  (1?))  constitute  a  oosvlete  solution  to  the  problem. 
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IT.  THB  VOS  ZEIPEL  PERTURB  ATIOH  TBCBHIQUB 

mm-tixm  sphdp  bcample 

to  far  the  equation  of  notion  of  a  mass  attracted  to  an  equi¬ 
librium  position  by  a  force  proportional  to  the  distance  has  been 
used  as  an  example  (i.e. ,  the  linear  spring).  Suppose  nov  that  a 
perturbation  tens  proportional  to  the  cube  of  the  displacement  is 
added  to  approximate  a  non-linear  spring  or  the  motion  of  a  pendulum. 
The  equation  of  the  system  is  nov 


••2  3 

mq^  +  k  <3x  +  ^<3x“° 


(23) 


e  q£ 


where  t  is  a  small  parameter.  It  is  necessary  to  add  the  term 

•e 

to  the  Hamiltonian  of  Eq.  (2)  so  that  the  canonical  equation  for 
becomes  equivalent  to  the  above.  That  is 


2  p  k 

B  Px  1  k2  2  Sc 

H  IT  2  k  Sc  +"TT 


<S*> 


so  that 

•  dH  .2  , 

Taking  the  results  of  section  II  In  expressing  the  new  Hamiltonian  I 


Px  •  *  §!  ■  m  m  *2  -  e  <1?  (®5) 


H(L,  X)  -  Hc  (L)  +  [h(L,  1)  -  Hq  (L)] 

vhere 

£  £ 

H-Ho  -  -Tf- 

The  portion  of  the  Hamiltonian  representing  the  perturbation  (H-Hq) 
contains  both  action  and  angle  variables.  The  unperturbed  portion 
contains  only  the  action  or  momentum  variable.  Using  the  definitions 
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of  L  and  X  and  the  unperturbed  solution  far  q^  of  lqs«  (15) 

g 

H(L,  t)  rn  kL  +  008^  X  (26) 

k 

The  first  step  of  von  Zelpel's  method  aonslsts  of  transforming  a 

Hamiltonian  H  (which  may  be  a  function  of  several  action  and  angle 

variables)  into  a  nev  Hamiltonian  H#  so  that  on*  or  more  of  the  an* 

# 

gular  variables  present  In  H  ana  eliminated  from  H  .  The  solution 

•  ** 

may  progress  in  stages  by  transforming  H  into  E  ,  etc . ,  until  no 
nngift  variables  axe  present  in  the  final  result .  Por  example,  in  the 
above  equation  H(L,  1)  would  be  transformed  to  -  ),  (the  dash 

indicating  the  absence  of  X').  i/  is  now  a  new  action  variable.  The 
resulting  H- J  canonical  equations  are 


V  m  ^5—  m  const.  i/  •  -  — —  •  0 
dl/ 

The  time  variation  in  i'  is  linear  (secular)  whereas  L#  is  a  constant 
of  the  motion.  It  may  appear  that  part  of  the  solution  has  been  lost, 
but  it  will  become  apparent  later  in  carrying  through  the  transformation 
from  H  to  H*  that  the  periodic  part  of  the  solution  is  obtained  in  the 
process.  The  final  results  of  the  problem  will  be  in  the  form 

X  ■  i'(t-tQ)  +  trigonometric  terms  periodic  in  X 

L  •  l'  ♦  trigonometric  terms  periodic  in  X 


cos  X  (using  X,  L  of  above) 


The  transformation  technique  is  described  as  follevs.  Assume  that 
H  »  +  (H-Hq)  may  be  expanded  in  an  infinite  series,  each  term 

having  a  coefficient  involving  some  power  of  the  small  parameter  l. 
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that  Is 

I  ■  ^  -♦  Bg  +  ••• 

share  the  subscript  refers  to  tbs  poster  ef  C  involved  la  the  ooe f- 

rt 

ficient  and  Hg(6  ) ,  tax  exwplo,  Is  designated  as  "of  order  2". 

Without  defining  its  aatheswtical  significance,  ^ 

t reduce  a  "deterulalng  Auction"  having  a  aev  variable  l'  and  the 
old  variable  i 

S  -  S(l/,  i)  -  so  +  Sj^  ♦  S2  +  *•• 


Cm  subscripts  have  the  some  meaning  as  before.  The  following  re¬ 
lationships  are  imposed  on  St 


L  m 


dS 

31 


3s 

te/ 


(27) 


It  is  desirable  that  l  •  i'uhen  6  •  0  (the  unperturbed  case) 

Bterefbre  is  arbitrarily  defined  as 
© 


S  «  L'  t 
© 


(28) 


The  canonical  transformation  is  to  be 


H(L,  *)  -  H*(l/,  -  ) 


or 


H0(L)  +  [l(L,  i)  -  Ho(L)]  -  H*  (L/)  4  <l')  4 

Jtaw,  using  the  relationships  from  Sqa.  (27)  and  (28) 

,  hsi  &s2 

H(L,  1)  *  H(L  +  £J-  4  gj~  +  *  •*,  1) 


(29) 


(30) 


•  Bc  4  %  +  Hg  + 
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and  expanding  the  Hamiltonian  of  Eq.  (30)  in  a  Taylor  aeries  about  1/ 
and  collecting  terns  of  corresponding  order  In  both  right  and  left 
members 


Zero  order 

VL 

')  -  <(L') 

<») 

First  order 

ds, 

~  +  Hjl/,  l)  -  h£(l') 

<32) 

Second  order 

K 

^  1  (hhY 

d!  2  bh'2  \bl/ 

(33) 

•S  &-*-> 

The  object  of  this  expansion  is  to  solve  the  resulting  partial 
differential  equations  (ordinary  differential  equations  in  this  spe¬ 
cial  case)  for  8^,  S^,  etc.,  depending  on  the  order  of  accuracy  de¬ 
sired.  Musn  this  is  done,  the  main  difficulty  of  the  problem  is  over¬ 
come  and  the  solution  follovs  at  once.  It  is  very  important  to  ob¬ 
serve  in  the  work  that  follows  that,  8^,  Bg,  etc.,  will  be  obtained 
by  solving  linear  equations.  That  is,  8^  vlll  be  determined  from 
Kq.  (32)  and  then  substituted  into  Xq.  (33)  which  la  linear  in  Sg. 
Similarly  the  third  order  equations  can  be  solved  from  an  equation 


linear  in  using  the  previously  determined  values  of  and  Sg. 
Consider  the  first-order  term 


H^L',  l) 


eos^  l 


(*) 


It  may  be  expressed  as  the  ran  of  two  terms,  one  constant  (l.e.. 
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"secular"  or  non-periodi c  in  X)  and  the  other  periodic  in  X.  These 
terns  are  easily  obtained  by  expressing 

OOS^  l  m  ^  +  £  CDS  2X  ♦  g  QOS  kl 


where 


H, 


or  altemativoly 


^(L',  »)  -  ♦  Bjj, 


-JSjr  *  fx  [”*  21  +  S  "•  kjt]  '  (M) 


_  &l/2  r*  _  k  .  3Rl' 

bi. •  7?  J„  “■  *“*i? 


n^p  *  ^  "  His 

The  first-order  Kq.  (32)  »ay  he  expressed  as 

&L  * 

—  ~  +  H.+  H,  ■  H, 

*/  ax  ^ 1  >  1 

and  separated  into  two  equations: 

/2 


“i.  •  <(L')  ’  ^r 


“ip 


&!.'  hi 


(*) 

(37) 


te  &H 

In  the  above.  — ®  stands  for  — -  with  L  replacing  L  after  the 

*L  . 

derivative  is  taken  or  in  our  socasplo,  einply: 

&L 


&/ 


*  k 
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Hence  from  Kq».  (35) 


^8,  or  l 

■  ""  a  rn  •  ( 008  21  +  008  bX) 

bt  acJ 


and 

(sin  21  ♦  gf  sin  kX)  ■*■  const.  (39)1 

i 

Hie  constant  of  integration  is  net  significant  to  the  problem 

since  we  are  only  concerned  with  the  derivatives  of  S^.  How  the 

part  of  the  solution  periodic  in  X  is  obtained  to  first  order 

(i.e.,  0(d)  fno«\»9*.  (27)). 

no  d(S^  ♦  S.)  -t/2  . 

L  ■  ■  11  . .  »  l/  «■  y*  (OOS  21  +  r  COS  4X)  (%0) 

ax  bt  at3  4 


i'  «  —  «  - fll  -  X  -  Q-  (sin  2/  +  J  sin  kl) 

-*,3  o 


L  m  i  +  ■  'k.  (sin  2X  +  i  sin  ht) 


If  first-order  accuracy  is  sufficient  in  calculating  the  periodic 
terns,  there  is  no  interest  in  determining  S^.  At  this  point,  the 
Hamiltonian  h*  is  fraa  Iqs.  (31)  and  (36) 

H*(l/)  -  H*  ♦  +  ... 


•  *  • 
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may  be  obtained  in  a  similar  manner  from  Eq.  (33) .  ft»l  8^, 
etc.,  vUl  all  be  periodic  la  4  cad  elnoe 

b2  H 


bL 


7 T 


*  ^si 

H^l')  ■  the  part  Independent  of  l  of -  - 

^  bh'  dX 


dSi 


Bov  — =  is  given  by  lq*  {38),  end  fm  Xq.  (fi6) 
bl 


&H1  2CL'  k  , 

-  ■  1  COS  £ 

bL'  k2 


so  it  is  necessary  to  determine  the  non-periodic  part  of 


2d2  L/3 


L  J  4  .  ,3  k 

008  X  {•£  -  008  4) 


considering  the  infinite  series  for  cos  X  (Bee  prior  to  Eq.  33) 


%COBk  l  -  ^ 


Isec 


8  # 
cos  £ 


“  "  fsl  +  B  +  T3s)  "  * 


sec 


Hence 


H, 


-s 


17  e2  L/3 


35 

m 


♦  ^ 

and  H  is  given  to  second  order,  o(£  ),  by 


H 


8k^  k5 


<**) 
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The  secular  part  of  the  solution  is  given  by  the  canonical  H-J 
equations 


dt  dl/ 


a 


/2 


dt  di 


<*3) 


Smnerlzlng  the  complete  perturbed,  solution  ve  have 


qJC(t)  ■  /\/ ^  cos  l  (kk) 

where 

l  m  t  (t-t^)  +  — — »  (sin  217  +  ^  sin  Ul*)  (hj) 

2kJ 


and 


i' 


i#(t-t0) 


-  rk  +| 


cl' 

7“ 


2  .  /2 


<t-tQ) 


L  *  L' 


(cos  2X* 


oos  ti') 


Hote  that  l  has  been  replaced  by  l'  in  Kqs.  (4o)  and  (fcl)  since  the 
difference  between  i  and  Ll  is  0(£).  The  coefficients  of  the  right- 
hand  terns  are  also  0(£),  making  the  discrepancy  small. 

Referring  back  to  the  original  variables  and  substituting  in 


L 


L' 


kSo2 

.t- 


{*7) 


Sco 


h 


(cos  2 l  + 


1 


cos 


w 


The  physical  interpretation  of  these  results  Is  as  follows:  <3^ 
may  he  Interpreted  as  the  actual  Initial  displacement  of  a  mass 
attached  to  a  nonlinear  spring  at  time  tQ.  The  displacement  after 
an  Interval  (t-tQ)  may  he  determined  from  the  simple  linear  spring 
equation  ty  considering  that  the  elements  (i.e.,  parameters)  defining 
the  linear  motion  are  time- varying  due  te  the  perturbation.  The  var¬ 
iation  of  the  linear  spring  parameter  defining  the  amplitude  Is 
given  by  gq.  (h8) .  The  variations  In  frequency  of  the  motion  are  in¬ 
herent  in  5q.  (h5). 

In  the  shove  exsaqjle,  the  perturbation  term  involves  only  the 
displacement  and  the  energy  of  the  system  is  conserved.  The  so¬ 
lution  of  gq.  (J*8)  shows  the  amplitude  of  the  oscillatory  motion  to 
have  no  secular  terms.  That  is,  it  does  not  Increase  or  decrease  with¬ 
out  bound  but  Is  represented  by  a  constant  plus  trigonometric  terms 
having  frequencies  which  are  multiples  of  the  fundamental  frequency 
of  the  system.  As  will  be  shown  later,  this  is  true  whenever  the  per¬ 
turbation  term  Is  independent  of  px> 

Bote  that  the  secular  part  of  the  solution  (gq.  %3)  is  accurate 
to  o(d  )  whereas  the  periodic  part  (gq.  kl)  is  accurate  to  o(&).  In 
theory,  the  solution  may  be  extended  Indefinitely  by  collecting  higher- 
order  terms  In  the  Taylor  series  expansion  of  H  and  solving  the  re¬ 
sulting  linear  partial  aquations  for  Sy  etc.  For  example,  con¬ 
tinuing  on  with  gqs.  (31)  through  (33),  the  third-order  terms  are: 


&/  hi 


bL'*  \bi 


bL'  bt 


-  H^(l') 


(*9) 


(Bote:  The  first  term  Involving  Is  periodic  in  l). 


has  been  determined  from  Sq.  (32)  and  Sg  may  be  determined  by 
finding  the  part  of  Bq.  (JJ)  periodic  la  i.  Shat  is,  the  periodic 
part  of 


gives  a  partial  differential  equation  which  may  be  solved  for  8g. 


<»> 


ItUUa  of  thla  caleulotloa  ox.  eorrlod  out  la  fifimUII  I  A. 

i.  ■ 

TOkCflcgrAJiyi  sigma 

Bo  tar  the  discussion  has  been  restricted  to  conservative  systeas 
where  there  is  neither  creation  nor  dlslpation  of  energy,  m  the  non¬ 
linear  spring  example,  the  perturbation  term  was  a  function  of  a  po¬ 
sition  coordinate  only  and  did  not  involve  the  time  explicitly.  Bext, 
let  us  consider  a  mechanical  system  to  be  perturbed  by  external  forces 
such  that  the  equations  describing  the  system  are: 
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vhere 

H  -  B(^,  px) 

Aa  In  tb»  preceding  oxmple,  let  as  stvpose  that  ««  vill  vlah  to 
transform  Kq*.  (52)  ao  that  they  art  la  terms  of  acaw  nsv  variables 
L,  X  aad  so  that 


and 


Kv  V  *  B'*L*  l) 


W  -4 

H  "  3T  +  x 


dL 

®-*ar*r 


(53) 


Tint  It  will  1m  necessary  to  express  X'  and  T/  la  terms  of  ths  asv 

variable*  L,  1.  Oils  transformation  of  variables  can  be  dona  la  the 

<*P,  dq_ 

following  manner.  First  multiply  Sqi.  (52)  by  aad  add. 


He  H 

ST  3T 


He  He 


ar  sr  • 


♦ 


If  sa  the  left  side  va  substitute 


He  He  dX  ■  H  dL 
dT  “  sr  e  +  sr  he 


H  He  di  A  He  di 

r -  3T3E  +  sr  i 


the  result  Is 


^ 

3T  3T  "  3T  1ST 


4j e  ai'  .r^iV^ 

H"3T  +  X3r +  Yar 


®*  teaekst  on  the  left  elds  is  a  Lagrangian  bracket  £  1,  L  J  and 
s^iale  unity  fbr  a  canonical  transformation  where  the  wint^n 
is  unchanged.^  (Bile  may  be  illustrated  by  carrying  oat  the  ebore 
indicated  operations  on  q*,  px  as  given  by  Bqs.  (15)  of  the  previous 
section) .  Bence 


—  **x  \ 

&-ar +  xar  4Tar 


di  _  ta'  .  _  ~'x 


Comparing  the  above  with  the  first  of  Kqs.  (53)  It  Is  that 


X^L,  t)  •  l{px,  q^)  ♦  T(px,  q^)  ^ 


m 


ip  ba 

If  Sqs.  (52)  are  aul tip  lied  by  jj-  ,  -  jy-  and  slnllar  operations 


carried  out  the  result  Is 


X  (L,  i)  ■  X  -rr-  ♦  I 


*>x  *Sc 


3T  *  x  3T 


(55) 


At  this  point  we  can  consider  an  exanple  of  a  non-conservative  sys- 
t«*  auefa  as  a  relaxation  oscillator,  let  this  system  be  described 
by  ran  der  Pol's  equation  which  is  of  the  fbza 


+  A 


Px  ■  0 
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Suppose  the  solution  to  he  periodic  of  the  fonn 


A*  V“  •<»-*.> 

px  -  ^  »  sin  »(t  -  to) 


(56) 


and  as  before  let 


‘-'14 


t  ■  »(t  -  tQ)  -  mt  +  B 


so  that  Sqs.  (56)  are  transfoxaed  to 


mr 

-  v  T  COB  1 


Px  ■  -  .y  2Lui  sin  1 


and  the  canonical  equations  are  nov 


•  _ _ o 

q*'sr 


(57) 


(58) 


(59) 


dl 

p  «  .  —2.  .  y 
x  d^ 


share 


FX  A  1  2  2 


X  -  0 


Upon  differentiating  Eqs.  (58)  it  is  clear  that 


b<± 

hi  *  m 


1 

3T  ~ 


Substituting  these  results  in  Eqs.  (5^)  and  (55)  results  in 

,  e(<^  -  i)px  r  2  1 

[  -  - . „;j - -  C|  tan  X  -  sin  X  cob  Xj 


OXJjj 


L 


,  £(<Jx  “  X)Px 

r  -  — — - -  -  e«B 


/  k  ^  2  \ 

2.  k  K, 

2 

(Sco  +  Sd) 

cos  X  -  cos  X  - 

V 

Referring  to  Eqs.  (53) ,  and  expressing  the  above  equations  as  a 
trigonometric  series,  ve  find  that 


dX 


+  t 


tan  X  - 


2 

V 

T 


sin  2X 


dl 


cos  21  - 


k 

So, 

nr 


cos 


(60) 


Bote  that  the  first  equation  vhlch  Involves  the  angular  rate  (and 
hence  the  periodicity)  contains  only  short  periodic  terms  whereas 


*7 


the  seoond  equation  involving  the  aaplitude  contains  a  secular  term 
pins  short  periodic  terns.  Taking  the  values  averaged  over  one  period 


1.  m  (t) 

Av 


LAv  "  T 


'  k 

So  2 

T"  "  So 


Using  these  average  values  in  order  to  obtain  a  first  approx¬ 
imation  to  the  solution,  ve  may  obtain 

So  “  2  Sd(1  "  -tt) 

(since  i  -  mq^  q^) 


Performing  the  integration 


(14) 


So 


So®  * e  (t-V 


<&) 


The  first  approximation  is  a  harmonic  oscillation 

S  "  So  0013 

where  the  amplitude  varies  according  to  Eq.  (6l)  and  the  frequency 
is  constant. 

The  preceding  example  vas  relatively  simple  because  of  the 
simple  form  of  the  Hamiltonian  Hq  (a  function  of  L  only).  Instead 
of  Eq.  53#  consider  the  case  where  the  canonical  equations  are 


where  I  »  Hq(L)  +  K^(1,L)  +  ♦•• 

the  subscripts  having  the  scone  significance  a a  in  the  example  of  the 
previous  section.  One  or  more  successive  canonical  transformations 
m«y  he  performed  by  the  von  Zeipel  method  so  aa  to  eliminate  the  angle 
variables  from  the  nev  Hamiltonian. 

Jbr  example 

H(L,/)  — »H#(l/  -) 

r  -  «.  x# 

&/ 

L'  rn  -  .  y' 

*l' 


Or,  if  two  transformations  are  involved 
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V.  APPLICATION  TO  dUCaKAL  MKBAHIOB 

The  non-linear  spring  and  Tan  der  Pol  equations  vers  chosen  be¬ 
came  of  ths  simplicity  involved  la  having  only  oas  tegree  of  free- 
don  whereas  perturbation  problem*  in  colsstlal  mechanics  usually  in¬ 
volve  either  three  or  flour  degrees  of  freedom.  As  la  the  preceding 
exanplea,  the  ultimate  objective  Is  to  express  the  effects  of  the 
perturbation  as  tine  variations  In  the  two- body  orbital  parameters. 

As  a  first  step,  the  Besiltoalsa  is  written  la  terns  of  a  convenient 
set  of  canonical  variables.  A  usual  procedure  flor  doing  this  is  to 
start  vlth  the  vs  11- known  and  basic  Delaunay  variables j  L,G,H  rep¬ 
resenting  momenta  and  l,g,h  representing  the  conjugate  angle  variables. 
Riot  le 

MOMEHTA  COOEDHAXKS 


G  -  Lv4-e8 


1  *  M  (mean  anoetaly) 
g  -  »  (long,  of  perigee) 


H  -  G  cos  1 


b  «  0  (node) 


vhers  a,  a,  1,  n,  m,  sad  T  are  the  familiar  non- canonical  Keplerian 
( 7 ) 

elements.' These  canonical  elements  may  then  be  regarded  as  a  key 
which  allows  the  Hamiltonian  to  be  written  at  once  and  leads  to  the 
canonical  H-J  equations.  A  very  important  point  is  that  the  Delaunay 
variables  may  be  transformed  in  a  prescribed  by  simple  algebraic 

manipulation  to  one  of  many  other  different  sets  which  may  be  more 
adaptable  to  the  particular  problem  (e.g. ,  Jbinca ri  variables  and 
variations  thereof) . ^ 
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In  celestial  Mechanics*  the  letter  7  ratter  than  E  la  usually 
used  to  denote  the  Samiltenian*  in  order  to  avoid  confusion  with 
Delaunay's  E. 

Then  letting 


will  give  the  correct  set  of  2n  -  6  canonical  H-J  differential  agna¬ 
tions  leading  to  the  unperturbed  two-body  solution,  since 


b¥ 

5T 


.  a? 

px  "ST 


when  P  ■  7q,  the  tine  derivatives  of  all  hut  one  of  the 
elements  are  zero*  i.e.. 


The  familiar  two-body  elliptical  equations  are  used  to  obtain  posi¬ 
tion  and  velocity  as  a  function  of  the  orbital  elements  and  time. 

In  writing  the  Hamiltonian  7  in  the  general  (i.e.,  perturbed) 
ease*  if  the  tins  la  present  explicitly  in  7,  an  additional  pair  of 
conjugate  variables  K,  k  may  be  introduced  for  convenience  as  a  de- 
•  vice  to  eliminate  t.^  The  solution  proceeds  by  successively  trans¬ 
forming  the  Hamiltonian  7  so  as  to  eliminate  the  angle  variables  t, 
g*  h,  or  k.  Each  transformation  may  involve  the  moral  of  one  (or 
more)  of  these  variables*  Such  a  transformation  might  go  as  follows 

♦Actually,  for  convenience  a  minus  sign  is  associated  with  7. 
that  is,  7  •  -  E  (the  Hamiltonian) . 
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F(L,G,H,  X,g,h)-~ f* (L#, O', 
using  S  ■  SQ  ♦ 

T  (Li,o',i',-,g,-)— f** {l',o#,hV,-,-) 

*  *  * 

using  S  »  80  +  ♦  ••• 

The  short  periodic  part  of  the  solution  involving  l  and  h  is  obtained 
by  solving  the  partial  differential  equations  for  S.  Similarly,  the 
long  periodic  terms  involving  g  case  from  S  .  The  secular  (linear 
time- varying)  changes  in  the  orbital  parameters  are  obtained  from 
the  resulting  H-J  equations. 


MLm, 

.  ^ 

dL*  < 

dt 

’  bL* 

dt 

hi 

.  ** 

m  v  ** 

do*  Sr 

dt 

’  do' 

dt 

Sg 

dh*  _  dF** 

as'  _  dr** 

dt 

W 

dt 

dh 

0 


TI.  THK  KHTLOTT  AMD  BOOCLIUBOrr  KSTHOD 


It  vu  shown  in  Section  III  that  the  effects  of  a  small  perturbing 
force  on  a  s  lapis  harmonic  action  could  be  described  or  approximated  by 
considering  the  elements  or  constants  of  the  linear  eolation  to  be  tlae 
variables.  Ibis  was  expressed  by  Bqs.  (22)  which  are  repeated  below 
using  the  more  cannon  notation  of  a  for  angular  frequency  instead  of  k. 

j£  -  T  ,ln  +  B) 

S  “  §  cos  (at  +  ») 

where 

X  -  X(px, 

The  Kcyloff-Bogolluboff  (K»B)  first-order  approximation  consists  of 

li  jm 

considering  ^  and  ^  constant  during  the  interval  frcn  zero  to 
2*/».  Accordingly,  the  objective  is  to  find  the  average  value  of  the 
functions  on  the  right-hand  aldea  of  the  shove  equations  during  this 
tlae  Interval.  It  stay  easily  be  shown  that  if  X  is  a  function  of  q^ 
only,  the  average  value  of  ^  will  be  zero  and  hence  A  will  be  a  con¬ 
stant. 

text,  consider  the  canonical  equations  as  discussed  in  Section  IV 


33 


where 

H  -  H(L,  Jt)  «  HQ  +  ^ 

t  *.  .  .  #  #  * 

H  -  H  (I/,  +  +  •••. 

h  -h. 


t  •  *(t-tQ)  •  «t  +  B 


dH 

3E 


-3T 


1  &H 


Referring  to  Eq.  (19)  f  It  la  clear  that  if  the  perturbation  tea 


X  is  a  function  of  only 


and  the  secular  or  constant  part  of  is  given  by 


*1*  *S 


c  *(<3*)  <1^  M 


(62) 


Using  the  t  -  o  solution 


SCO 


cos  jt 


3^ 


Equation  (62)  lacawi 


2  x  a -i 

6q*o  X^S»  °°*  ^  ila  *  *J*  «*  (63) 

di' 

whets  ilii  dtsssy  variable  of  integration.  therefore,  Jg-  la  given 
to  first  order  in  d  by 

<*> 


It  la  interesting  to  ccnpane  the  above  procedure  with  that  of  the 
K-B  first-order  approximation.  First  It  should  be  acted  that  by  def¬ 
inition 


di 

n 


-  0)  ♦ 


dB 

d£ 


Using  a  prime,  as  in  the  preceding,  to  indicate  the  perturbed  value, 
the  K*B  first-order  approximation  is  determined  by  averaging  ^  over 
the  period  ~  so  that 

V 


lr "  ® +  ms^)  "h  J*  x^v> 008  1)  cos  i  di 


(65) 


It  is  dear  that  for  the  special  case  where  XCq^)  is  of  the  form 


X(qx)  -  e^°  -  liq^  oos  if 


where  a  is  same  constant 
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©sen  o f  Xq.  (6k)  beooees 


sad  Sqb.  (6k)  and  (65)  give  th •  result 


(«) 


0 


Tbs  preceding  may  be  extended  for  tboee  cum  where  xfq^)  may  be  rep¬ 
resented  as  a  power  series  or  polynomial  In  q^. 

10X-0QHSXK7ATI7X  8Y8TB18  OSIKO  THE  K-B  METHOD 

In  Section  17,  the  Van  der  Pol  equation  Is  considered,  using  a 
canonical  system  of  equations.  The  ease  solution  nay  be  obtained 
using  the  K-B  first  approximation  by  assuming,  as  before,  that  the 
solution  Is  of  the  form 


^  -  V  cos  *(t-to) 

kxm  m  clxo  *  sla  «(t-to) 

Then 

X(V  p^)  -  -  (i  -  q^)  -  (l  -  cos2  /)  »  sin  1 

"  Sco  *  [C1  “  TT  }  Bin  1  “  TT  8in  31] 
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Applying  the  general  relationships  of  Eqs.  (22)  and  finding  the 
Minn  valus  orer  ths  interval  i  ■  »(t~to)  -  2* 


dX 

3C 


«  »  + 


£ 

Tn 


•In  X  cos  X 


2 

■^2  sin  3X  cos 


1 


di 


0 


sin 2X 


2 

«XD 

TT 


■In  3X  «in  X 


di 


These  rasults  agree  with  the  average  or  secular  terns  as  derived  fxm 
Sqs .  (6o) . 
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TU.  OMTCLPSIOHS 

Sh  precedir^;  Material  primarily  concern*  the  use  of  canonical 
variable*  and  the  Bamllton-Jacobi  equations  in  finding  approximate 
solutions  to  non-linear  mechanics  problems.  Sine*  von  Zoipal'*  per¬ 
turbation  technique  la  similar  in  many  respect*  to  the  Kcyloff- 
Bogoliuboff  method,  some  discussion  and  coaparieon  of  the  two  is 
given,  particularly  since  the  latter  i*  popular  fbr  solving  problems 
in  the  fields  of  engineering  and  physics  whereas  the  fhrmer  is  popular 
in  astronomy.  Both  methods  allow  the  solution  of  a  complicated  non¬ 
linear  problem  to  be  separated  into  a  preconceived  desirable  fhrm 
(e.g. ,  separation  of  periodic  terms  into  convenient  groups). 

As  pointed  out  in  Ref.  (7)>  the  use  of  a  canonical  syetsm  of 
equations  penults  a  transformation  to  be  applied  to  any  number  of 
terms  in  the  disturbing  function  simultaneously,  sad  the  Hamiltonian 
is  obtained  directly  in  the  course  of  finding  the  determining  function. 
Moreover,  the  new  canonical  variables  resulting  from  the  transformations 
are  related  to  the  original  non- canonical  set  by  simple  formulas  arising 
from  the  unperturbed  solution.  Ibe  solution  may  he  carried  out  theo¬ 
retically  to  any  desired  degree  of  accuracy  in  terms  of  powers  of  a 
small  parameter  6.  Them  a  user's  standpoint  it  is  convenient  because 
the  solution  is  obtained  in  a  very  methodical  way  involving  no  great 
mathematical  difficulty. 

Another  interesting  point  la  that  the  method  of  solution  does 
not  contribute  erroneous  secular  tome  to  the  solution.  As  pointed 
out  by  Poincar^^  and  others,  the  classical  perturbation  technique 


3S 


(dm  to  Foiaeon)  giv»«  rim  to  limar  tlm- varying  taras  vtxlch  do 
not  hav*  physical  mailing  (a.g. ,  a  Umar  tlm-lnoraaalag  aapUtuds 
In  tho  imlimar  spring  psoblan)  * 
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Appendix  a 

HiagguQBZacR  aoLunoi  qr  pow-ldbar  spruc  prqbum 

As  outlined  in  Sect  leu  IV,  the  solution  nay  be  carried  out  to 
higher  orders  of  accuracy  by  evaluating  successively  the  partial 
differential  equations  resulting  from  the  Taylor  Series  exp  ana  Ion 
ef  the  Hamiltonian.  8^  has  boon  detoxmiaod  from  Sq*  (32)  as 

gr  <2 

ff  «  -  —-nr-  (sin  21  +  1/8  sin  kt) 

A  lar3 

Sg  may  bo  determined  by  finding  the  part  of  Sq.  (33)  periodic  in  1. 
That  is,  the  periodic  part  of 


Hirst  we  accumulate  the  following  relationships 

to/ 

dV 

_ lx  m  0 

to/2 


to/ 


cos 


(from  Sq.  3*) 


3 r 


1 1  1  (cos  1  -  3/8)  (from  Sq.  38) 

k'3 
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Hance,  free  the  above 

— 5  -  +  periodic  part 
*X  k° 


[eos8  X  -  £  eos*  ij 


the  trigDametrie  identity  fbr  ooe^  X  is  given  prior  to  ®q.  (35)*  the 

fi 

identity  for  cos  X  is 

COS8  X  -  yjjfl  [cos  8X  +  8  008  6x  ♦  28  008  4X  +  56  cos  2X  +  35J 
fnoa  which  we  nay  derive  the  fallowing: 

^  [00s  2X  +  yg  oos  4X  ♦  ^  cos  6X  +  cos  8xJ 


fV3 


[sin  2X  ♦  ||  sin  4X  ♦  ^  sin  61  ♦  ^ 


2-  /2 

•  'g"  [sin  2X  ♦  ^  sin  ti  +  ■<jg  sin  6x  +  ysi 

4k 


ds,  8s 

iote  that  ^  etc.,  being  constrained  to  be  purely  periodic 

*  * 

in  X,  cannot  contribute  to  or  Hg  respectively.  ft>  determine 
,  we  refer  to  Sq..  (49)  which  is  repeated  belov: 


»tm3 


hi 


part  independent  of  t  of 


1 

I 


sr*r 


m  part  Independent  of  i 


_  t  h 

of  — K  OOB 

k 


h 

OOB 


Carrying  out  the  indicated  multiplication  and  retaining  only  the 
oonstaat  or  non-periodic  terns,  we  obtain 


Dae  secular  solution  of  Iqs.  (h3)  nay  now  be  modified  to  Include 
terns  of  O(C^) 


Doe  periodic  portion  of  the  solution  (Sq*.  ho  and  hi)  mgr  be 
similarly  modified  to  include  terns  of  0(£^).  Dae  man  angular  rate 
is  expressed  in  terns  of  the  initial  displacement  by  using  the 
substitution 
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